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ABSTRACT 

Many quantum computation and communication schemes require, or would significantly benefit from, true sources 
of single photon on-demand (SPOD). Unfortunately, such sources do not exist. It is becoming increasingly clear 
that coupling photons out of a SPOD source will be a limiting factor in many SPOD implementations. In 
particular, coupling these source outputs into optical fibers (usually single mode fibers) is often the preferred 
method for handling this light. We investigate the practical limits to this coupling as relates to parametric 
downconversion, an important starting point for many SPOD schemes. We also explored whether it is possible 
to optimize the engineering of the downconversion sources to improve on this coupling. We present our latest 
■ results in this area. 

1. INTRODUCTION 

' The advent of photon-based quantum cryptography, communication, and computation scheme j- 1 l 2 l 3 l 4 l 5 l B l 7 l 8 l IJ l lu l 

has increased the need for light sources that produce individual photons.^ An ideal single-photon source would 
produce completely characterized single photons on demand. Since all of the currently available sources fall 
far short of this ideal, many new schemes have been proposed for creating single photon on demand sources 
(SPOD). Many of these schemed- 1 ^ 13 * rely on optical parametric downconversion (PDC), because it produces 
photons two at a time, allowing one photon to herald the existence of the other. In a previous work we proposed 
one such scheme for a multiplexed single photon on-demand source that increases the probability of single photon 
emission while suppressing the probability of multi-photon generation.^ Most PDC based schemes (including 
ours) require that the PDC output be collected into a single spatial mode defined by an optical fiber. In order for 
these PDC schemes to produce a reliable SPOD source, it is essential that the optical collection system efficiently 
gathers and detects the herald photon and sends its twin to the output path of the system. 

Ch ' 

Various theoretical models have been d eveloped to predict how the collection efficiency of PDC light in a 
"two-photon single mode" can be improved ! 14 * 15 ^ 16 ! and in many cases these models have been experimentally 
proven to give a better coupling efficiency. In particular, it has been shown that the position and size of the 
pump beam focus affect the shape of the PDC outpulf^and hence the coupling efficiency of the PDC into a given 
spatial mode. In addition, a more detailed 

worlPI 

recently showed how crystal length and walk-off also affect 
the coupling efficiency for a pulsed broadband pump. (However, in that recent work the attempts to increase the 
single mode fiber coupling efficiency were centered on having a thin crystal and a tightly focused pump beam, 
which is achieved at a cost of overall source brightness.) The work done to date on coupling efficiency gives 
some guidance for increasing coupling efficiency. However, there is still not a practical "recipe" to quantify and 
maximize the coupling of PDC light into single mode fibers. 

In this paper we present a method for analytically calculating the coupling efficiency of PDC light into single 
mode fiber. Specifically, we calculate the spatial and angular overlap between type I PDC output and the two 
single modes defined by optical fibers. The collection efficiency is evaluated in terms of pump waist, collecting 
mode, crystal length, and the imaging conditions of the optical collection system. To ease the computational 
difficulties associated with this problem, we make the approximation of near-perfect phase matching. We quantify 
the limitations introduced by this approximation and then show by example how we can obtain better coupling 
efficiency by properly shaping the pump beam and collection modes for a given crystal length. Although our 
approximations limit the length of crystals that can be analyzed using this method, this formalism does allow us 
to study crystals with lengths beyond what is considered the thin crystal limit. As a next step in this research 
we plan to use a more exact model with a numerical evaluation that allows us to model the efficiency over a 
wider range of parameter values (expecially longer crystal cases). 
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2. THEORY 

To determine the coupling efficiency between the parametric downconversion output and two single collection 
modes (defined in our setup by single-mode optical fibers) we must calculate the spatial overlap of the two-photon 
field of the PDC source with the modes selected by the fibers. The first step in this procedure is to calculate the 
two-photon fielcEl given by 

a 12 = (oi4 +) 4 +) i^), a) 

where |0) is the vacuum state and e[ + ^ and E^ are the positive- frequency portions of the electric field evaluated 
at positions ri and Y2- The two-photon state at the output surface of a PDC crystal oriented with its face 
perpendicular to the z-axis is given by 



J d 3 k s J d 3 h J d\ p J d 2 p £ dzE^K^^+^ll^l^) 



(2) 



where S is the cross sectional area of the crystal illuminated by the pump, p is the transverse distance from 
the z-axis, and L is the length of the crystal. We designate the component of a wavevector that is parallel to 
the crystal face by k, so that k., = k Z jZ + Kj. The subscripts j = s, i, p indicate the signal, idler, and pump. 
The quantities Ak z = k Pz — k sz — k\ z and Ak = referred to as the longitudinal and transverse 

wave- vector mismatch, respectively. The pump beam's angular spectrum defines its transverse distribution via 
the Fourier transform 

1 f ,2 



E P (P) = ^J d z Kp E p (K p y K * >>. (3) 

This analysis assumes that the pump propagates with negligible diffraction effect inside the crystal, so that E p (p) 
is independent of z. 

The electric fields in Eq. are evaluated outside the crystal, while the wavevectors in Eq. J2J) are evaluated 
inside the crystal. To evaluate Eq. it is convenient to write Ak z in terms of the angular frequency, u>, and 
the index of refraction, n(u>), since the z components of the wavevectors are discontinuous at the crystal surface. 
Thus, the z component of a wavevector is 

/ / nY/< jV, i \ 2 

(4) 

where c is the speed of light. We assume that the pump beam has a narrow angular spectrum (transverse 
wavevector distribution) and that the signal and idler are observed only at points close to the z axis. In 
addition, we assume that only a narrow range of frequencies are collected because narrow bandpass filters are 
placed the signal and idler optical paths. The central frequencies defined be these filters are specified by fi s 
and Qi. (In our notation we denote central frequencies chosen as system parameters by capital f2 and frequency 
variables by lower case ui.) We then expand Ak z around these central frequencies to obtain^ 

Ak z « Ak° z +Duj- D'uj 2 + 2ns(a)acos0i ^ + 2^(^)^008 02^ " 2 n p (9 p ,n p )n p K P' (5) 




where we have 
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In obtaining these expressions we have assumed type-I phasematching with perfect frequency matching, i.e. 
O p = Oj + f2 s where fl p is the pump angular frequency. The angle between the pump propagation direction and 



Figure 1. Schematic of collection setup. 



the crystal optic axis is specified by 9 p , and the PDC photons emission angles are given by 9\ and 9 2 . We have 
developed the angular frequency as u s = Q s + lo and uji = Oj + u>' , with u' = —u>. Note that in the degenerate 
case we have D = 0. 

We assume the PDC light is collected into two single-mode fibers using identical lenses (with focal lengths /) 
both placed at a distances d from the fiber tips and Md from the crystal face, where M is the magnification. (We 
assume that the fibers define single Gaussian modes. See Fig. 1 for a schematic of this setup.) In this situation 
we can write the biphoton field as 

A 12 ( Pl ,p 2 ) oc J diJ J d2 Ks J d 2 Ki J c\ 2 k p E p (k p ) J d 2 ^ dz (7) 

x e^ z+AK -^H i (K i ,p 1 )H s (K s ,p 2 )T 1 (Lo)T 2 (Lo), 

where -ffi(/«i, pi) and H s (k s , p 2 ) are the transfer functions of the lenses in the optical path of the idler and signal 
and p\ and p 2 are the transverse coordinates measured from the z\ and z 2 axes, respectively (see Fig. 1), at 
the imaging plane of the lenses. The functions t\(u>) and t 2 (uj) define the spectral transmittances of interference 
filters placed in the two collection paths. To simplify the form of the transfer functions, we use the paraxial 
approximation and consider an ideal lens with infinite aperture. Under these assumptions the transfer functions 
are given by 

Hi(K U pi) oc e ^.-Pi el c e d 2 K I 2 /(2n I ) 

H S (K S , P2 ) OC e *MK s -p 2 Jced^l/(2n a ) 

where e = 1/d + 1/ (Md) — 1/ f provides for non-perfect imaging caused by precision limits in the setting of the 
lens- fiber- crystal distances. 

The pump field transverse distribution is defined to be Gaussian with a waist of W p , and we assume that 
the transverse crystal size is large relative to the pump beam, so we can take the cross section S to be infinite. 
We also assume perfect transverse phase matching, i.e. k p = k s + k^. These assumptions allow the k p and p 
integrals in Eq. Q to be trivially performed 



J s d2 p(J d 2 *p E p (K p )e^ P^ e -«(".+«0-P oc E p {k, 
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We next assume that the longitudinal mismatch is small (Ak z L -C 1) and that we have near-perfect imaging 
(ced 2 K 2 i /(2fl Sy i) <C 1)- In this case, we can expand the exponentials in the Eq. (JJJ and use the result in Eq. J§J 
to obtain 



(10) 



l + ^Afc 2 i-i(Afc 2 L) 2 



The spectral filters defined by t± and t 2 are centered around the wavelengths f2 s and f2; are both assumed to have 
a rectangular bandwidth B. We limit the expansions to first order and assume degenerate type-I phase- matching 
(i.e. D = 0). The integrals yield 



Ai 2 {px,p 2 ) oc BA° 2 (p 1 ,p 2 ) 
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where the zero order solution for the biphoton field is 
A 12 {px,p 2 ) oc j^Ep 



and the coefficients are 



5(pi-p 2 ) / dwTi(w)T 2 (w). 
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With Eq. (|ll|l in hand we are prepared to calculate the coupling of the biphoton field into the single mode 
fibers in our setup. To accomplish this we write A 12 as a coherent superposition of guided modes in the fiber as 
suggested in Ref. ED 

A 12 ( Pl , P2 )= A 12 lW ' lm l pt, mt (p 1 M m (P2) (14) 

I'm' ,1m 

The collection efficiency can then be calculated by 
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where C\ 2 , C\, and C 2 measure the square of the overlap between the biphoton field and the collection modes, 
and can be used to calculate singles and coincidence counting rates 
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j l d2p2 Al2 ( Pl > P^lmiPl) 
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We assume the guided mode is a gaussian denned by 



rio(p)=y!^«p 



(17) 



where Wf is the width of the collection fiber. The waist of the collected mode imaged back to the crystal is 
given by W, where W = MW/. The collection efficiency is then obtained by using Ea. pi|l and (jl5f) 



V12 ^ V12 



o 1 + Wx\ 
1 + \L 2 X2 



where the collection efficiency at the zero order of the perturbation theory (very thin crystals) is 



^12=4 
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and the corrections are 



Xi = 



X2 



Jd 2 p 2 Jd 2 p 1 A° 12 (p 1: p 2 )v* 1Q (pi)p* 10 (p2)(a + Pai + 7a 2 - fa 3 ) 
Jd 2 p2 / dV^i 2 (Pi 5 P2Mo(PiM (p 2 ) 

/ d 2 |Q 2 |J d 2 ^! A° 12 {pi, p 2 )(fl (pi)(a + /?ai + 7 a 2 - ^Qa)j 
Jd 2 p 2 |/d 2 piA° 2 (p 1 ,p 2 )^* (p 1 )| 2 



(18) 



(19) 



(20) 



We can perform these integrals using the expressions developed above to obtain 
2v 8(3 -2^ + 87 



Xi 



X2 = a 



W 2 W 2 + 2 W2 
2 , 12^ 2 , 16(/3 + 7 )^ 2 1 



1 



1 



W 4 



W 2 W 2 



W 4 (W 2 + W 2 ) 2 



„, ) + 4a( ^ + (/3 + 7)( >^ + >v 2 + w 2 



(21) 
))■ 



3. EXAMPLES 

To illustrate the use of Eq. (|18|) we consider a Lithium Iodate crystal with L = 2 mm. The crystal is pumped 
with 458 nm light with 9 p = 36° and the collection modes are arranged to collect degenerate PDC emission at 
Q\ = 6 2 = 1.5° with a bandwidth of B = 2 x 10 13 s _1 (about 10 nm). The magnification is fixed at M — 20 and 
we assume that the collection mode waist occurs at the crystal with perfect imaging (e = 0). Figure 2a plots the 
collection efficiency versus the pump waist for different collection mode waists. Figure 2b plots the term Ak z L, 
and Figs. 2c and 2d plot the correction factors Xi-^ 2 /4 and X2^ 2 /4, respectively. Notice that the correction 
factors become significant for small pump waists (i.e. broad angular spectra of the pump). 

Figure 3 plots the collection efficiency versus the pump waist for four different crystal lengths: L = 0.1 mm, 
2 mm, 3.5 mm and 4.5 mm, at various mode waists (all are in the valid range of the approximation). Note that 
for the case of L = 0.1 mm (which satisfies the thin crystal approximation) the optimum collection efficiency 
can be obtained for all the collection modes by picking a sufficiently large pump waist. However for the longer 
crystals this is not the case. For example, if the collection mode is defined with W — 20 /im in the case of 
L = 4.5 mm then the maximum collection efficiency that can be obtained is 7712 ~ 0.9 regardless of how the 
pump waist is varied. However if we choose a collection mode with W = 50 /xm the collection efficiency is 
r)i2 ~ 1 when Wp > 100 /im. (We limit the analysis here to crystal length of 4.5 mm because longer crystal push 
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Figure 2. (a)Plot of the collection efficiency 7712 calculated versus the pump waist for various collection waists: W = 
50 fim, 60 fim, 70 /im, 80 /J,m, and fOO /im. The crystal length is L — 2 mm. (b) Plot of the longitudinal mismatch Ak z L 
in the same condition as case (a), (c) Plot of the correction factor x?£ 2 /4. (d) Plot of X2L 2 /4 — the lines are so close 
that the individual collection mode waists are not labelled. 




Figure 3. Plot of the collection efficiency 7712 calculated versus the pump waist for various crystal length from left to 
right L =0.1, L—2, L=3.5 and L=4.5 mm. Each graphics is evaluated for collection mode waists W = 20 (jm, 25 /im, 
30 /im, 35 /im, 50 /im. 

the validity of the approximations made earlier.) From Fig. 3 we conclude that we can take advantage of the 
higher signal levels from long crystals without sacrificing coupling efficiency provided that the pump waist and 
collection mode are scaled to appropriate values. 

In Fig. 4a we plot the collection efficiency versus the crystal length for different magnifications (M=15, 20, 
25, 30, 40) while maintaining a constant pump waist (W p = 50 /im) and fiber size (W/ = 1.5 /im). Figure 4b is 
the same as Fig. 4a but in the non-perfect imaging condition e = 0.4m _1 , which reduces the collection efficiency. 
Notice that as the crystal gets longer larger collection modes give a better coupling efficiency. Even though a 
longer crystal length reduces the collection efficiency obtainable for a given pump size, these longer crystals may 
still be desirable since they produce more overall signal. In addition, if the pump beam size is also be adjusted 
(as was done in the Fig. 2) it is still possible to get good coupling efficiency with a longer crystal. 

To verify that Eq. i|18|l is valid for the parameters that we have chosen, we must consider the numerical values 
of the corrections and verify that our approximation is in its range of validity. This is shown in Fig. 5 where the 
correction factors ~x\L 2 /A and X2L 2 /A are plotted versus the crystal length for different magnifications (M=15, 
20, 25, 30, 40), maintaining the pump waist at W p =50 and assuming the fiber core of W/=1.5 /im. It is 
clear that for the chosen parameters for crystals longer than about 8 mm the corrections become non-negligible. 



4. CONCLUSIONS 

In the interest of optimizing the brightness of two-photon sources for use in emerging applications of quantum 
information, we have presented an analytical calculation of the collection efficiency of PDC light into single mode 
fibers. Our calculation was performed for the case of long pump pulse and type I PDC. These experimental 
conditions are similar to those used in many SPOD sources. We demonstrated that for crystals that are longer 
than those required by the thin crystal approximation (but still within the approximation made in the model) 
the maximum collection efficiency can be obtained at larger pump waists. However, as longer crystals are used 



h> = 22.5 |im 




L (mm) 



Figure 4. (a)Plot of the collection efficiency 7712 calculated versus the crystal length at W p =50 /im for various collection 
mode waists. Heavy lines are calculated for perfect imaging (e = 0) and narrow lines are calculated for the case of 
e = 0.4m" 1 . 
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Figure 5. Plot of the correction factors Xi^ 2 /4 and X2L 2 /4 versus the crystal length for different collection mode waists, 
maintaining the pump waist at W p =50 (jm. 



the maximum achievable collection efficiency decreases. The approximations made in obtaining the expression 
for coupling efficiency limit the range of parameters that can be studied using this method. It is particularly 
important to pay close attention to the validity range of the approximations to assure meaningful results, as this 
has not been emphasized in some other works. Future work will move beyond the simple two-photon efficiency 
analysis to include an overall signal estimate that accounts for the effect of finite available pump on the total 
signal received. Experimental tests will be conducted to verify these analytical predictions. 
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